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ABSTRACT 

The inner crust of mature neutron stars, where an elastic lattice of neutron-rich nuclei 
coexists with a neutron superfluid, impacts on a range of astrophysical phenomena. 
The presence of the superfluid is key to our understanding of pulsar glitches, and is 
expected to affect the thermal conductivity and hence the evolution of the surface 
temperature. The coupling between crust and superfluid must also be accounted for 
in studies of neutron star dynamics, discussions of global oscillations and associated 
instabilities. In this paper we develop Lagrangian perturbation theory for this problem, 
paying attention to key issues like superfluid entrainment, potential vortex pinning, 
dissipative mutual friction and the star's magnetic field. We also discuss the nature 
of the core-crust interface. The results provide a theoretical foundation for a range of 
interesting astrophysical applications. 



MOTIVATION 



The crust of a mature neutron star shields the high-density 
fluid core from the thin atmosphere. It provides a heat- 
blanket that determines the link between the neutrino- 
driven cooUng of the core and the evolution of the observed 
surface temperature (Gudmundsson 1983). The physics of 
the crust is also of key importance for neutron star dynam- 
ics. The presence of a superfluid component, and the asso- 
ciated rotational vortices, is central to our understanding 
of pulsar glitches (Espinoza et al 2011). The coupling be- 
tween the neutron-rich nuclei in the crust and the superfluid 
neutrons also affect global oscillations involving the crust, 
as in the case of the quasi-periodic oscillations observed 
in the tails of magnetar flares (Andersson, Glampedakis & 
Samuelsson 2009). On the one hand, our understanding of 
the crust physics is quite good. The equation of state for 
matter has been modelled in detail. In particular, we have 
a clear idea of how the composition changes with density 
which allows us to work out the elastic properties of the lat- 
tice of nuclei (Chamel & Haensel 2009). This also provides 
us with the fraction of neutrons that have dripped out of 
nuclei and which may considered free to move relative to 
the lattice. However, our understanding of the dynamics of 
the coupled crust-superfluid system is quite poor. This is 
somewhat surprising given the relevance of the involved is- 
sues (e.g. in the context of glitches), but it is nevertheless 
the case. 

To make progress we need to develop a theoretical 
framework that allows us to model the dynamics of the 
inner-crust region. Such a model has to consider the key 



physics, like elasticity, the magnetic field and superfluidity. 
It must account for the presence of superfluid vortices, po- 
tential pinning and mutual friction. It is also important that 
the model is adaptable, in order to faciliate further develop- 
ment concerning, for example, finite temperature effects and 
heat propagation. This is obviously a major challenge. The 
aim of the present work is to develop a "complete" model for 
linear perturbations of a (moderately) realistic neutron star 
crust. Working within the Lagrangian perturbation frame- 
work (the natural way to view the perturbed crust, and a 
necessity if one wants to consider the stability of the system 
from the formal point of view (Friedman & Schutz 1978)), 
we provide a two-component model that accounts for the 
key physics. The model also complements previous work on 
the outer-core region (Andersson, Comer & Grosart 2004) 
by discussing the role of the superfiuid entrainment and the 
magnetic field (Glampedakis & Andersson 2007). We also 
consider the conditions that need to be implemented at the 
crust-core transition. The theory is developed to the point 
where the results can be applied to a range of interesting as- 
trophysics problems. Having said that, there is massive room 
for future refinements as our understanding of the relevant 
physics improves. 



2 THE TWO-FLUID MODEL 

The conditions in the inner crust, with an elastic lattice of 
increasingly neutron-rich nuclei immersed in a neutron su- 
perfluid has (at least) two dynamical degrees of freedom; 
the superfluid neutrons may flow relative to the lattice. In 
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essence, this is an example of a two- "fluid" system, although 
in this case one of the components is also affected by elas- 
tic restoring forces. Our model for this system is based on 
the variational approach to multi-fluid dynamics developed 
by Prix (2004) and Andersson & Comer (2006) (see also 
Haskell, Andersson & Comer (2011), which represents the 
Newtonian limit of the fully relativistic convective varia^ 
tional model designed by Carter (1987), see also Andersson 
& Comer (2007). The elastic sector builds on the relativis- 
tic model developed by Karlovini & Samuelsson (2003), and 
the inclusion of superfluidity follows the strategy set out by 
Carter & Samuelsson (2006); Carter, Chachoua & Chamel 
(2006). 

2.1 Hydrodynamics, magnetic fields and elasticity 

In order to model the conditions in the inner crust of a 
neutron star, we take as the starting point the equations for 
two-fluid hydrodynamics (Prix 2004; Andersson & Comer 
2006; Haskell, Andersson & Comer 2011). Assuming that 
the individual components are conserved (and working in 
a coordinate basis where vectors are represented by their 
components, with indices i,j,k as usual), we then have the 
usual conservation laws for the number densities rix, where 
X is a constituent index labeling the components. 



dtn^ -h Vi(nxvi) = 



(1) 



Note that a repeated species index x docs not imply sum- 
mation, in contrast to the spatial indices like i for which the 
Einstein summation convention applies. In the outer core of 
a neutron star, the distinction between the two components 
is fairly clear. On the one hand, we have the supcrfluid neu- 
trons. On the other hand, we have a conglomerate of (most 
likely superconducting) protons and electrons. On scales 
larger than the electron screening length and timescales 
longer than the inverse of the plasma frequency (McndcU 
1991; Glampedakis, Andersson & Samuelsson 2011), these 
form a single, charge neutral, fluid. The two-fluid model sim- 
ply distinguishes the neutrons (represented by rin, say) from 
the protons /electrons (given by rip). The situation in the 
crust is similar yet different. At densities beyond neutron 
drip, some neutrons remain bound in nuclei but there is 
also a "gas" of free neutrons. The assignation of neutrons to 
each component follows from the equation of state, once the 
nature of the ions in the lattice is established (Chamel & 
Haensel 2009). However, if we consider the dynamics of the 
system it is not clear to what extent the "confined" neutrons 
are able to move (Carter, Chamel & Haensel 2006). This 
depends on how strongly bound they arc, to what extent 
they can tunnel through the relevant interaction potentials 
etcetera. The upshot is that one can choose to work in differ- 
ent chemical "gauges" . This issue has been discussed in de- 
tail by Carter, Chamel & Haensel (2006) (and we will return 
to it later). The choice of chemical gauge affects the inter- 
pretation of the involved quantities (number densities, etc), 
but the two-fluid model remains unchanged conceptually. To 
make a distinction from the outer core problem, we will refer 
to the two components in the crust as "free" neutrons, with 
density rif, and "confined" baryons, represented by ric- This 
notation may represent a slight bias towards to description 
advocated by Carter, Chamel & Haensel (2006) and Carter 
& Samuelsson (2006), but at this point we basically want 



to keep the options open by not linking the discussion too 
much to established results for the outer core. 

The "free" neutrons can flow relative to the crust lattice 
once the system cools below the transition to ^So neutron 
superfluidity (see Andersson, Comer & Glampedakis (2005) 
for typical transition temperatures) . We then have two cou- 
pled equations of momentum balance. Assuming that the 
large scale system comprises a sufficient number of quan- 
tized vortices that macroscopic averaging is meaningful (this 
should, indeed, be the case for all astrophysical systems of 
interest) the two momentum equations can be written (Prix 
2004) 

(2) 

where x = {f, c}, the velocities are wj, the relative velocity 
is defined as wl^y = vl^ — Vy and jl^ — /ix/rrix represents 
the chemical potential scaled to the nucleon mass (we will 
assume that the neutron and proton masses are equal, rric = 
rrin = me). The mass densities are given by px = msrix and 
$ represents the gravitational potential, which means that 
we have 



: 47rG(pf -I- pc) 



(3) 



and, finally, the parameter Ex encodes the entrainmont. The 
entrainment can be expressed in terms of a single parameter 
a such that (Prix 2004) 

PxSx = 2a (4) 

To close the system, wo need to provide an equation of state. 
In the present formalism, the equation of state takes the 
form of an energy functional E, the functional form of which 
determines the chemical potentials 



dE 

and the entrainment parameter 

dE 



(5) 



(6) 



The forces on the right-hand side of (2) can be used 
to represent various other interactions, including dissipa- 

tivc terms (Andersson & Comer 2006; Haskell, Andersson 
& Comer 2011). If we focus on the conditions in the neu- 
tron star crust, then we need to account for elasticity and 

the large-scale magnetic field. The latter is described by the 
usual electromagnetic Lorcntz force. That is, we have 

/r = Kju-pB" (7) 



Eliminating the total current with the help of Ampere's law, 
i.e. J' = (c/47r)e'^'°VjSfc, this becomes 



ft 



(8) 



In order to use this in (2) we also need to know ric, i.e., to 
what extent the baryons are affected by the magnetic field. 
It is natural to assume that all baryons that are confined to 
nuclei are involved, but the exact meaning of this is (as we 
will discuss later) somewhat fuzzy in a dynamical situation 
in regions of the stars where some of the neutrons are free. 

The elastic force is different in that it only serves to 
restore deviations from a relaxed state of the lattice. This 
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means that it is natural to discuss elasticity at the linear 
perturbation level. We will do this later, in Section 4.2. In 
that discussion, we will assume that the background con- 
figuration is relaxed, in which case there is no leading order 
elastic force in (2). This assumption is not quite realistic; the 
crust of an astrophysical neutron star is likely to be strained 
due to the regular electromagnetic spin-down of the system. 
It is, in principle, straightforward to account for this strain 
but in the interests of clarity we have chosen not to do so 
here. 

Before we move on, it is worth noting that the incor- 
poration of elasticity requires us to track given "fluid" ele- 
ments relative to the relaxed configuration. This motivates 
us to use a Lagrangian framework. This approach is, of 
course, advantageous for a number of reasons. In particular, 
if we are interested in considered rotating neutron stars. Our 
model for the inner crust builds on the two-fluid perturba- 
tion framework developed by Andersson, Comer & Grosart 
(2004), adds the magnetic field according to the analysis of 
Glampedakis & Andersson (2007) and provides a model for 
elasticity which represents the Newtonian limit of the theory 
developed by Carter & Samuelsson (2006). 

2.2 Superfluidity and vortex dynamics 

Let us turn to our attention to the superfiuid aspects of 
the problem. Doing this, we note that (2) accounts for the 
presence of a (macroscopically averaged) vortex array. In or- 
der to discuss issues concerning, for example, vortex pinning 
and mutual friction it is useful to consider (2) in more detail. 
Following Glampedakis, Andersson & Samuelsson (2011) we 
introduce the momentum 

p* = m (^ul + ctwf^ (9) 

In a superfluid, the momentum arises as the gradient of the 
condensate wavcfunction. The upshot of this is that the su- 
perfluid is irrotational. However, this is only true on the 
microscopic scale. On the scale of hydrodynamics, the su- 
perfluid can rotate by forming vortices and when these are 
averaged over the system mimics a system with bulk rotation 
(as evidenced by (2)). The rotation is, however, quantized 
and the vorticity is given by 

^ leiifeVjpl = nv4 (10) 
rn 

where riv is the number of vortices per unit surface area and 
(with a unit vector along the direction of the 
vortex array and k = h/2m « 2x 10~^ cm^ s~^ the quantum 
of circulation). 

Prom equation (10) one can derive the equation that 
governs the vorticity. Assuming that the vortex density is 
conserved we have 

9tnv-|-Vi(nv4) = (11) 

where vl, in fact, defines the macroscopically averaged vor- 
tex velocity. The fact that the vortices move with vl also 
means that (in terms of the Lie derivative Cv^ along vl) we 
have 

{dt + jO^,) K = (12) 
Given these relations, it follows that 

dtWl + eijuV' (e'^'-^Wlvl,) = (13) 



In order to make contact with the macroscopic descrip- 
tion, we now rewrite the relevant Euler equation, c.f. (2), 

as 

dtp\ + Vi [jif - yt)f -I- Pj- j - meijkVfWi = ff/rif (14) 

This expression makes it clear that, in the absence of vor- 
tices and external forces, the superfiuid motion follows from 
the gradient of a scalar potential. Moreover, it is now easy to 
compare (13) and (11). As discussed by Glampedakis, An- 
dersson & Samuelsson (2011), we then find that the models 
arc consistent provided we account for the "Magmis force" 
on the right hand side of equation (14). This force takes the 
form 

f( ■ k k 

— = nvtijkiiHvt - i)v) (15) 

Pf 

and an equal and opposite force will act on the vortex array. 
In the simple case of a single condensate at zero tempera^ 
ture, force balance on the vortices requires them to flow with 
vl (here and in the following wc ignore the inertia of the vor- 
tices (McndcU 1991)). In a more general situation, wc can 
still use the above strategy to account the forces that act on 
the vortices. We "simply" solve the force balance equation 
for the vortices for vl. and use the result in (15). In the case 
of resistive scattering off of the vortex cores, e.g. by phonons, 
this leads to the usual representation of the vortox-mcdiatcd 
mutual friction (Andersson, Sidery & Comer 2006) . We will 
now consider this problem in the context of the crust. 



3 MUTUAL FRICTION AND VORTEX 
PINNING 

To discuss the various vortex forces, we takes as our starting 
point the equation of force balance for a single vortex; 

t'^''ki{vi-vi) + n{vi-vi) + r = Q (16) 

This accounts for (i) the Magnus force, (ii) a resistive friction 
Eissociated with the normal component (e.g. nuclei, electrons 
and phonons in the crust) , with coefficient TZ, and (iii) a gen- 
eral force which we leave unspecified at this point. This force 
will later be taken to represent the "pinning" of vortices to 
the nuclei in the crust. Note that the different terms all have 
the dimension of velocity. In order to obtain an expression 
for the force per unit length of a vortex, as required in the 
macroscopic Euler equations, we need to multiply by ptKriv 
This step assumes that the vortices form a recti-linear array 
(the usual Abrikosov lattice). It is the simplest set-up, essen- 
tially since it make the averaging over vortices trivial, but 
there is no guarantee that this situation prevails in a real 
system. In particular, it may be relevant to worry about the 
formation of vortex tangles and superfluid turbulence (An- 
dersson, Sidcry & Comer 2007; Pcralta et al 2006). We will 
not consider this problem here. Neither will wc account for 
contributions like the vortex tension or the elasticity of the 
vortex array. These effects are readily incorporated in our 
framework (see Haskell (2011) for a recent discussion of the 
elasticity of the vortex lattice) but we leave thcin out in the 
interest of clarity. Finally, we assume that there are no force 
contributions along k', the direction of the vortex axis. 

Solving for the vortex velocity in the standard way, we 
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find 



where 



r = e'^'=K,^r+^' (18) 

Now, the reaction force acting on the normal component 
will be 



J c T-j / i i 



(i + 7e2) 

l + 'R? ' 



(19) 

where we have defined the projection orthogonal to the vor- 
tices; 

±}= 5] - k'kj (20) 

The force acting on the neutrons will naturally be equal and 
opposite; 



(21) 



For later convenience, it is natural to introduce basis vectors 
along the macroscopic relative flow (in the plane orthogonal 
to the vortex). That is, we use 

where 



and w 
together with the decomposition 



= ail + a±e''''' KjWk 



(22) 
(23) 
(24) 



This means that we get 



(1 +7^-' 



= [TZ^w + a\\TZ — Ox] e'^'^kjWk 

- [Tlw + a±n + ay] (25) 

Expressing the result in the usual form (Andersson, Sidery 
& Comer 2006) we have 

A 

PfK 



(26) 



where 



and 



Beff 



I + 11"^ 



9 a\\Ti. — ax 



(27) 



(28) 



The first term in each bracket represents the standard mu- 
tual friction. The second terms illustrate how a pinning force 
may be accounted for in the macroscopic multi-fluids model. 
It is worth noting that this model can also be applied to the 
problem of (potentially strong) interaction between neutron 
vortices and proton fluxtubes in the outer core of a neutron 
star, c.f. Glampedakis, Andersson & Samuelsson (2011). 



3.1 Perfect pinning 

Having discussed the general model, we are equipped to con- 
sider the limiting case of perfect "pinning" . The interaction 
between the vortex lines and the crustal nuclei may be strong 
enough to "pin" the vortices and force them to move with 
along with the crust (Donati & Pizzochero 1978). This has 
profound implications for the macroscopic dynamics of the 
system. Given that the vortex lines are no longer free to 
move, the superfluid neutrons cannot spin down (or up). 
Hence, a lag wifl build up between the two components as 
the crust slows down due to magnetic braking. When this 
lag develops, the Magnus force will tend to push the vortices 
out (or in), c.f. (15). Eventually, the force will be strong 
enough to overcome the pinning and break the vortices free. 
This leads to a transfer of angular momentum, that could 
explain large pulsar glitches (sec Sidery, Passamonti & An- 
dersson (2010) for a recent discussion). Vortex pinning may 
also have a severe effects on neutron star precession. By 
acting as a gyroscope, the pinned vortices are expected to 
lead to extremely short period precession, of the order of 
the rotation period (rather than the several months to years 
period expected from a typical crustal deformation) (Jones 
& Andersson 2001; Link 2009). While this general picture 
is supported by a range of theoretical models, we are still 
quite far from a detailed understanding of the nature and 
strength of vortex pinning. However, for the present study 
it is sufficient to assume that a pinning force is acting. 

Let us assume that there is "perfect" pinning, v\ = Ve- 
in this case the equation of force balance for a single vortex 
takes the form 



- «fc) 



(29) 



As long as the system remains below the unpinning limit, 
this provides us with the required pinning force J-^ . Given 
this, we flnd that the force acting on the neutrons is: 



fi ijk ^ / c f 



(30) 



while the reaction force on the charged component will be 



PfK 



J c ijk ^ / c 



fx 
■ Vk) 



(31) 



The main conclusion from this exercise is that, if the lag 
between the two components is sufficiently small to allow 
us to consider the vortices as perfectly "pinned" , the exact 
form of the pinning force does not appear explicitly in the 
equations of motion. 



3.2 Vortex creep 

Let us now consider the situation where the lag between 
the two components is close to the critical value for unpin- 
ning. This regime is of great physical interest as it is at the 
heart of many theoretical models for pulsar glitches. If one 
assumes that parts of the system are always slightly subcriti- 
cal, one may allow for a population of thermally excited vor- 
tices to unpin randomly and transfer angular momentum to 
the crust. This is usually refered to as "vortex creep" (Alpar 
et al 1984). Describing this behaviour is clearly a challeng- 
ing task, both from the microscopic and the macroscopic 
point of view. On the one hand, there have been efforts 
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to calculate the pinning "force" and the barrier that the 
thermally excited vortices would have to overcome to unpin 
(Link 2009). On the other hand, there have been attempts 
to incorporate the concept of vortex creep in a macroscopic 
hydrodynamical description, by assuming that only a frac- 
tion of the vortices, on average, participates in the dynamics 
(Jahan-Miri 2006). 

Here wc adopt a phenomenological approach aimed at 
exploring the hydrodynamics of the creep regime. We start 
by noting that vortex creep would correspond to motion such 
that vl « v*. Assuming that (24) describes the "pinning 
force" completely, i.e. that there is no component along k', 
we can rewrite (17) as 

Vcr =vl-vl = Y^p^ [ail''?- - Wcf - a±] wlf 

+ ^ ^j^2 ["II + (™<:f + ^] f'^'^'^J ^fc' (32) 

From this expression we see that there are two ways of en- 
forcing vortex creep. One would be to let 7?. ^ oo. This 
model has recently been used in studies of precession and the 
unstable r- modes (Glampedakis & Andcrsson 2009). How- 
ever, with the "pinning" force explicitly in the problem we 
have another option. Focussing on the 72. <C 1 case, which 
should be relevant in the crust (Feibelman 1979; Bildsten & 
Epstein 1989), we can demand that 

a\\R. — Wci —a± = v\\ < Vc (33) 

and 

a\\ + (wcf + a±) 7? s < Wc (34) 
With these definitions the creep velocity is given by 

Ucr = V||W' + V±e'^'' KjWk (35) 

and we have 

Vcr =V\\+v'i (36) 

This model is, obviously, more complicated since we now 
have three coefficients to specify; we need TZ, ay and a±. 
However, this provides more flexibility and could allow us 
to, for example, consider a specific form for the "pinning" 
force or, indeed, the creep rate Vct- An important point is 
that, in this model you do not have to have strong drag, 
7?, ^ 1, to effect pinning. This may be particularly relevant 
if wc want to make our neutron star precession models more 
realistic,see Link (2003). 

Let us conclude by writing down the force that enters 
in the hydrodynamics. After some straightforward algebra 
we have 

a|[ ~ v±+ U||7?. (37) 

and 

a± fs —w — V|| + v±TZ (38) 
This means that the force becomes 

—2- fti -v±w^ + {w + v\\)t^^kjWk (39) 

which means that 

fc ~ v]_ + [w + v\\)^ K vo^ (40) 



The required hydrodynamical force follows once we multiply 
by ptUvK. 

It is important to remember that the "pinning" force 
discussed here does not describe the realistic interaction be- 
tween the vortices and the nuclei. We have discussed a purely 
phenomenological model which allows us to describe the mo- 
tion when vortices are not yet completely free so that the 
drag force caiinot be approximated as linear in the veloci- 
ties (Link 2009). Our simple "pinning" model allows us to 
consider the dynamical implications of this regime. 



4 LAGRANGIAN PERTURBATIONS 

So far, we have discussed the general conditions that prevail 
in the inner neutron star crust and a two-fluid "hydrody- 
namics" model that accounts for the key features. We will 
now take an important step towards astrophysical applica- 
tions by developing a framework for Lagrangian perturba- 
tions of this system. The aim is to provide a model that 
can be applied to a range of important problems in neutron 
star dynamics, from pulsar glitches to the gravitational-wave 
driven instability of the r-modes and magnetar oscillations. 
These problems are all naturally approached within pertur- 
bation theory. Moreover, they all represent scenarios where 
the sensitive interplay between the crust, the superfluid and 
the magnetic field is expected to be important. 



4.1 The unentrained two-fluid problem revisited 

Since the two-fluid problem has two dynamical degrees of 
freedom it is natural to introduce two distinct Lagrangian 

displacement vectors ^* (Andersson, Comer & Grosart 
2004). In order to distinguish between these displacements, 
we use variations Ax such that (for any, scalar or vectorial, 
quantity Q) 

A^Q^SQ + C^^Q, (41) 

where 5 represents an Eulerian perturbation. The perturbed 
continuity equations, c.f. (1), then take the form (Friedman 
& Schutz 1978) 

AxTix = -rixVi^i — >■ 5nx = -V,(nx^x) ■ (42) 

This means that the equation that describes the perturbed 
gravitational potential is 

V^5$ = 47rmBG(5nx + Suy) = -A-irmsGV i{n,,ii + riy^) 

(43) 

Considering the simplest case of vanishing entrainment 
and no "external" forces, // = fi = 0, the results of Ander- 
sson, Comer & Grosart (2004) show that 

{dt + r„J Axuf + Vi ^Ax$ + AxAx - ^^xWx^ = (44) 

After some algebra, this leads to 

+ $xViVj$ - + ViA^p,^ = . (45) 

Here, the Lagrangian perturbation of the chemical potential 
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can be written (with y 7^ x) 

(46) 

using the fact that /ix = fi^{ni, ric). Hence, we arrive at the 
following final form for the perturbed Euler equations; 

dUt + 2viVjdt& + (viVjfG + Vi5$ + eiVi V,- ($ + Ax) 



9nv 



V,(ny^J) 



(47) 



Andersson, Comer & Grosart (2004) demonstrated how 
one can proceed further and derive useful conserved quan- 
tities, extending the single-fluid analysis of Friedman & 
Schutz (1978) to the two-fluid arena. The importance of the 
obtained canonical energy Ec stems from the fact that it can 
be used to assess the stability of the system. In order for the 
evolution to be dynamically unstable, i.e. for a perturbation 
to blow up in absence of additional forces, we must have 
Ec = 0. A secular (viscosity or radiation driven) instabil- 
ity requires Ec < 0, provided that the energy lost through 
dissipation is positive (which makes sense). A particularly 
nice feature of the analysis of Andersson, Comer & Grosart 
(2004) was the proof that the standard instability criterion 
for gravitational-wave instabilities, that a normal mode of 
oscillation becomes unstable when its pattern speed changes 
sign (eg. when an originally backwards retrograde mode in 
a rotating star becomes prograde (Andersson 2004)), holds 
also in the two-fluid problem. This had previously been as- 
sumed to be the case, but there was no formal proof. In 
the following, we will not attempt to address the issue of 
the canonical energy for more complex systems; our focus is 
entirely on the perturbed equations of motion. A stability 
analysis for these equations would be interesting, but as this 
may well be prohibitively complicated we leave the problem 
for future considerations. 



4.2 Accounting for elasticity and the magnetic 
field 

The equation that represents that crust dynamics must ac- 
count for both elastic and magnetic contributions. The for- 
mer leads to the well-known contribution 



Ac (/r/Pc) = -V^aij 

pc 

where the shear tensor is given by 



(48) 



(49) 



(here one should not confuse the shear modulus jl with the 
chemical potentials /Ux). It is important to keep in mind that 
these expressions are only valid for unstrained background 
configurations. 

In the case of the magnetic field, we need the Lagrangian 
perturbation of (8). The required results have already been 



derived by Glampedakis & Andersson (2007), so we simply 
restate them here. The perturbations are determined from 



which leads to 



and 



Ac ( — 1 = 

.Pc 



AcB' = -B'V^d 



(50) 



(51) 



Ac Si = BjVi^ - BiiViii) + B'Vjii (52) 
Finally, using 

Ac(V,BO = Vj(AcBi) - BtVjVi^l (53) 
we obtain from (8) 

Ac (/r/Pc) = -^[V,(AcBO - Vi(AeB,)] (54) 
47rpc 

It is straightforward to express the magnetic perturbations 
in terms of the displacement but since expressions that 
result are quite involved we have decided not to do so here. 

4.3 Entrainment 

The elastic and magnetic forces take relatively simple forms. 
However, the fluid part of the problem becomes much more 
complex if we consider the generic situation where the en- 
trainment is not vanishing. In the unentrained case the two 
equations of motion are coupled chemically through the 
equation of state (e.g. through various interactions) and 
gravitationally since variations in the number density of one 
fluid affect the gravitational potential and hence the motion 
of the other fluid. In contrast, the entrainment parameter a 
encodes how the internal energy of the system depends on 
the relative velocity of the two fluids. This usually leads to 
a stronger coupling of the components. Since the entrain- 
ment encodes the effective dynamical mass of each matter 
constituent, it will affect most scenarios that involve the su- 
perfluid regions of the star. 

Including the entrainment, the Euler equations take the 
form (2). However, since 

viVj (exwf'') = i^^wl^) - e^wY'^ivi (55) 
these can be rewritten as 



{dt + £„J « + ExwD + Vi + jSx - = 



(56) 



We want to consider Lagrangian perturbations of this sys- 
tem. The derivation follows the same approach as in the un- 
entrained case, described by Andersson, Comer & Grosart 
(2004) and summarised in Section 4.1. 

First we note that the continuity equations and the 
Poisson equation are not affected by entrainment, so we can 
still use equations (42) and (43). Secondly, perturbing the 
Euler equations we have 

(^t -h J [Axwf + Ax (ex^D] 

+ Vi f Ax* + AxAx - = (57) 



This follows immediately since the Lagrangian variation 



Lagrangian perturbations of the neutron star crust 7 



commutes with dt + Cv^- Most of the terms in this equa- 
tion were already considered in the unentrained problem. A 
key difference now is that J/ix depends on the entrainment. 
This is obvious since the chemical potential, /ix, is the par- 
tial derivative of the energy functional, E, with respect to 
the number density, rix. Since E depends on the entrain- 
ment the Eulerian variation of the chemical potential also 
depends on the entrainment. In general, we have 



5Ax 



9/ix 



where 



and 



givmg 



W-B V(^nx/„„,„2 TUB 



6w = 2w^^5wy^ 



(60) 



5Ax = 



(nxCi) 



^) V, (n,^^:) + ^^xwf 5< (61) 
At this point it is worth noting that 

- (^^ - ii) ^jvi + <v,e; - c^v,t.;x 

= {dt+C,^)(fy-Q)~C^,A; (62) 

The only other piece of equation (57) that was not 
present in the unentrained problem can be written 

(ft-h£„JAx (ex^r) 

= £x {dt + ) Axwf + {9t + ) AxSx (63) 
The first term follows easily from 



(64) 



The second term in (63) is worth discussing in more detail. 
We first consider 

AxEx = &x + dVj£x (65) 

Using the definition for Ex, and the perturbed continuity 
equation, we find 

2 

Px 

In the general case, the entrainment parameter a is a func- 
tion of the two number densities, e.g. nt and ric, and w^. 
This means that 



AxEx = — \5a + Vi faCi 
Ox L V 



(66) 



where we have defined 



da 



(67) 
(68) 



We can also use 

V^a = A^V^nf + AcViUc + 2AwwfV^wl( (69) 
This means that we can write Sa as, 

5a = -A{V^ {ntiPj - A^V^ [n^C) 
+2A^wf [dtii - dtii + ^'cVi^^ - vlViii - CNivi + Qv4 



(7 



After some algebra, we finally find that 
2 f 

Ax£x = —{ (a - ^xTlx) ViQ, - AyUyViiy 

+2^^«;f [dt - - vtVi^ + v;v4 - (4 - &) ViV^] 

Wc can now combine the above results to get a general 
expression for the right-hand side of (63) in terms of the 
two displacements. However, this expression will be rather 
lengthy and may not be particularly useful. In most situa- 
tions of interest a reduced version should suffice. In princi- 
ple, one may consider different simplifying assumptions. The 
most drastic would be to consider the entrainment param- 
eter to be uniform. The natural way to achieve this would 
be to take a = constant^. Then we have A^ = Aw = 0, and 
the equations simplify greatly. In fact, from (71) we are only 
left with 



2a 

AxEx = Vi^x = SxVi^x 

Px 



(72) 



A more realistic model would be based on an expan- 
sion for small relative velocities (Comer & Joynt 2003). One 
would expect Wyx to be small in most cases, so it makes 
sense to use the approximate equation of state 

S(nf,nc, w^) w Eo{n!,nc) + Ei{n!,nc)w^ (73) 

In this case, we simply have 

a = Ei (74) 

and it is obviously the case that Aw = 0. Hence, we have 

AxSx = — \{a- ^xTix) Vi^i - A'^y Vi^l (75) 

Px L -I 



This expression completes our analysis of the perturbations 
of the inviscid problem for the coupled crust-superfluid sys- 
tem. 



4.4 Perturbing the mutual friction 

The various contributions associated with the supcrfluid vor- 
tices add further complexity to the problem. As an illustra- 
tion of this we will focus on the mutual friction, which follows 
from equation (26). To complete the perturbation equations, 

we need 



Ax I ^ 



-Ax = (^B'efft^'^nvKjWk — KUvBeSw'' 



(76) 



^ It is worth noting that £x = constant is only consistent for a 
uniform density model. 
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where it is worth recalling that we defined Wf = tIvK' in 
Section 2.2. 

Prom equations (27) and (28) wc see that the pertur- 
bations of the coefficients Beff and B'^g will (in general) re- 
quire a knowledge of TZ, ay and a±. The perturbations of 
these quantities can, obviously, be treated in the same way 
as a in the previous section. However, in this case any sim- 
plifications would rely on an understanding of the detailed 
microphysics. In addition to these quantities we need the 
variation of the magnitude of the relative velocity, w. Deter- 
mining this quantity is straightforward. 

The main new piece of information required for the mu- 
tual friction is the perturbed vorticity. In general, when the 
vortices arc not moving with either of the macroscopic flu- 
ids, we need to perturb (10), after solving for the vortex 
velocity as in Section 3. The procedure is relatively straight- 
forward, but as the final expressions are messy, and not very 
instructive, we will not work out the details hero. Instead we 
consider the two extremes of free and pinned vortices. 

In the first case, when the vortices are free so that = 
Uf , we see that (10) leads to 

{dt + C^i) W't + WiiVy^) = . (77) 

Perturbing this, it is quite easy to show that (in the case of 
a stationary and axisymmetric background) 

[dt + C^) [AiWi + WKV.CD] = (78) 

We need the trivial solution to this equation, which means 
that we have 

AfWl = -WliVj^i) (79) 

Pinally, we perturb (11) (which is completely analogous 
to the continuity equation for n^) to get 

Af-riv = -rivVjCf • (80) 

Given these results, and the discussion in Section 3.1, 
it is easy to work out what happens when the vortices are 
(perfectly) pinned. In that case, we have = vl and as a 
result we find that 

AcWf" = -wK^di) (81) 

and 

AcTiv = -n^Vjii . (82) 



5 THE CRUST-CORE INTERFACE 

In order for the developed perturbation framework to be 
useful for neutron star astrophysics, we need to consider 
the crust-core interface. This region is known to be impor- 
tant for a range of problems, especially since the associated 
viscous boundary layer may provide efficient dissipation of 
large-scale flows in the core. We will not consider the vis- 
cous problem here, but it is worth keeping in mind that a 
key issue concerns to what extent the velocity perturbations 
are continuous across the interface. If they are not, then vis- 
cosity works to smooth out the discontinuities (over some 
relatively short length scale) leading to damping of the bulk 
motion. The magnetic field may play a similar role. When a 
magnetic field penetrates the interface, discontinuities would 



induce Alfvcn waves which would effect an efficient coupling 
between the crust and the core. 

Another issue arising in this context is the potential ap- 
pearance of nuclear pasta, that the nucleons form non-trivial 
topological clusters (e.g. rods or plates) rather than spherical 
nuclei arranged in a Coulomb lattice. Accounting for these 
structures is, in principle, straightforward once the proper- 
ties of the various phases are understood. In most cases one 
would expect the system to remain "isotropic" on macro- 
scopic scales owing to the fact that the pasta structures will 
"freeze" in a random fashion on some smaller scale, and the 
"fluid model" arises from a larger scale average. However, 
there are cases where these structures may be aligned on 
macroscopic scales (for instance due to a strong magnetic 
field or the existence of an ordered array of vortices). Then 
we may need to consider non-isotropic elasticity. The com- 
putation of the microscopic input parameters (equation of 
state, entrainment, shear modulus, et c.) is very complicated 
in the pasta phase. It may, for instance, be that the densities 
exhibit discontinuous jumps across the crust-core interface, 
which coxiki lead to discontinuities in the velocities [see e.g. 
equation (97) below] and thus to enhanced viscous damping 
as discussed above. These arc very important issues, but for 
simplicity we will ignore them in the following, assuming an 
isotropic solid and continuous densities across the crust-core 
interface. 



5.1 Chemical gauge 

As before, we consider a system where a charge-carrying 
component is coupled to a neutral superfiuid. Furthermore, 
we assume that the superfiuid extends across the interface. 
At this point we have to return to the issue of the "chemical 
gauge" . That is, we have to discuss the physical meaning of 
rii and ric- In developing the model, we have taken the view 
that rtf represents the neutrons that are not confined to nu- 
clei in the crust, while ric represents all protons as well as the 
confined neutrons (making up the ions in the lattice). This 
view leads to a natural description of the elastic and mag- 
netic forces. However, it does not lead to a straightforward 
connection to the core, where one would usually distinguish 
all the neutrons, rin, from the protons, rip. Problems arise 
from the fact that the analysis requires variables that are 
"meaningful" across the crust-core interface. In principle, 
the problem would be more straightforward if wc were to 
use a two-fiuid model based on rin and rip also in the crust. 
The downside to this would be that we would then have to 
reconsider the Lorentz force and the elasticity contributions. 
After all, some of the neutrons will be associated with the 
nuclei and hence should be affected by the crust motion. 

Focussing on the generic case, we will connect the stan- 
dard two-fluid model for the core with the crust model we 
have developed. This forces us to consider the relevance of 
the chemical gauge and serves to clarify some of the key 
issues. 

The chemical gauge choice relates to the neutrons that 
are considered "free". The issue is subtle since, in a dynamic 
situation, even the neutrons that are associated with the 
nuclei may be able to tunnel through the relevant interac- 
tion potential. This maies concepts like the atomic number 
somewhat hazy. In general, one may introduce a new basis 
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such that 



n't = nl, + {1 - ac)np 



(83) 



where Oc (which we will take to be constant in the following, 
a good approximation at the level of the individual fluid 
elements) accounts for the fact that some of the neutrons 
move with the (crust) protons. We also have 



(84) 



Given these relations, it is easy to show that the neutron 
momentum is independent of the chemical gauge (Carter, 
Chamel & Haensel 2006) . This follows immediately from the 
definition of the momentum ( Andersson & Comer 2006) ; 



Pi 



dC 

dni 



(85) 



where C represents the relevant Lagrangian. Hence, we have 



Pi =Pi 



It also follows that 



/in 



/if 



(86) 



(87) 



However, these results also show that in general we must 
have Vf ^ Vn and ef ^ e^. Finally, in order to consider the 
vortices across the interface, it is natural to assume that 

Wj = Wi (88) 

The behaviour of the vortices may, of course, be more com- 
plicated than this but it makes sense to first consider the 
simplest "reasonable" model. 



5.2 The background configuration 

Following Glampcdakis & Andersson (2007) wo represent 
the moving interface by a level set of a scalar function / 
which can be extended in a smooth fashion. Expecting the 
interface to move with the charged component, we require 



[a+£„J/ = , 
from which it is easy to show that 

[dt + Cv,]V^f = 



(89) 



(90) 



In other words, the gradient Vi/ is constant in the frame 
moving with v^. Prom this, it follows that the perturbation 
Ac/ satisfies 



[dt + £„JAc/ = 



(91) 



The trivial solution to this equation is Ac/ = 0, which es- 
sentially means that a fluid element at the original surface 

remains at the perturbed surface. The normal to the surface 
can obviously be taken to be Ni = Vif, and hence we have 
the unit normal 

Ni = Vif/N , where N = |V/| = [g'^VifViff^ (92) 

This means that 

[dt + C^,]Ni = Ni [dt + C]N . (93) 

which shows that, even though Ni is not preserved by the 
flow, any change in the unit normal is parallel to the normal 
itself. We will use this fact later. 

These considerations are quite general. However, in the 



problem of interest we may restrict ourselves to configura- 
tions (at the unperturbed level) that are stationary and ax- 
isymmetric. These assumptions mean that we have dtNi = 
and Nivl^ = 0. The latter represents a no-penetration condi- 
tion, simply stating that (in the background configuration) 
the core fluids do not migrate into the crust. 

In order to obtain the interface conditions, we identify a 
small cylinder of fluid aligned with the normal to the inter- 
face, Ni. Integrating the various equations over this small 
volume, we will be able to deduce the relevant conditions 
to impose. Carrying out this exercise we need to make sure 
that the equations we consider are valid in both the crust 
and the core. Given this, it is natural to take as our starting 
point the conservation of baryon number. Assuming that the 
problem is stationary, the core equation 



Vi (rij^vl + Upvl^ = 



matches to the crust result; 



: 



(94) 



(95) 



Integrating these over the small volume, we see that we 
should impose 



Ni (unvl, + Upvlj = Ni (niVf + ricvl^ 



(96) 



at the interface. Noting that there are no chemical gauge 
issues concerning the protons, the corresponding conservar 
tion law leads to (assuming that the densities are continuous 
across the interface); 



Ni (vi - vj) = 



(97) 



Given this, and the fact that the total number density is 
given by n = rip -h rin = nf -|- nc, we can rewrite the first 
condition as 



finiVi (vl, - = n{Ni (vf - vl^ 



(98) 



The final conditions (97) and (98) are, of course, trivially 
satisfied for an axisyrnrnetric system since NiV^. = 0. 

Moving on to the momentum equations, it is natural to 
work in the frame moving with the protons/crust. After all, 
the volume that we integrate over is fixed in this frame. We 
can also safely treat any relative flow as planar, since the 
volume we consider is arbitrarily small. These assumptions 
simplify the analysis greatly. 

Let us first introduce the total momentum flux; 



TTi = P^Vi + PyVl 



(99) 



where (x, y) is either (n, p) or (f, c), depending of whether 
we consider the core or the crust. Combining the Euler equa- 
tions in the appropriate way we find that (Andersson & 
Comer 2006) 

dt-m + V,- {vi'K^ + vl^nl) + Vip + pVi$ = V'Tij , (100) 
where the pressure p is defined such that 

Vip = n^ViH^ - aViW^y . (101) 

X 

The elastic and magnetic stresses are accounted for in Tij . 

Let us consider the problem on the crust side of the 
interface (the core follows simply from letting f — >^ n and 
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c — >■ p). In the frame moving with the crust we have = 
(obviously), and we also need to replace Vf Wfc- This leads 
to 

7r| ^ pf (1 - , (102) 

and the momentum equation takes the form 

dtm + Vj [pf(l - EfXwf + 51p\ +pV,$ = V'Tij , (103) 

We now integrate this equation over the small volume strad- 
dling the interface. As long as the total density is continu- 
ous across the interface, the gravitational potential and its 
derivative will be smooth and therefore the corresponding 
integral vanishes as we let the volume shrink. This exercise 
tells us that there will be no local force associated with the 
interface as long as 

L J crust 

= Nj U(i - s^)wi^wr + sip - Ti] (104) 

L J core 

This condition is quite general. In particular, it needs to 
hold also on the perturbativc level. As far as the background 
configuration is concerned, we obviously have Njv^^ = 0, 
which means that we have (representing the change in a 
given quantity across the interface by (. . .) (Glampedakis & 
Andersson 2007)); 

Nj(Sip-T,^)=0 , (105) 

These are the usual traction conditions. 
For the vertical component see that 

{p) = N^N^T,j) , (106) 

while the horizontal components lead to 

jj' {Tij) = , (107) 

where we have defined the projection (orthogonal to the nor- 
mal) 

_l"= (?" - iV'iVV (108) 

It is relevant to note that the pressure may now be affected 
by the presence of a relative flow. 

If we assume that the background configuration is such 
that the crust is relaxed, we only need to account for the 
magnetic stresses. Then we have 

Tij^-gr,^ + ^B.Bj . (109) 

We also know that ViS' = 0, which implies that we must 
have 

{NiB*} = . (110) 
In this case, condition (106) leads to 

(P+^>=0, (111) 

while (107) becomes 

(^NjB'^ _L" {B,)^0 . (112) 

Combined with (110), this shows that if the magnetic field 
penetrates the interface then all components of the back- 
ground field Bi must be continuous. 



To complete the analysis of the interface, we need one 
more condition. We obtain this condition from the momen- 
tum equation for the neutrons. This choice is natural since 
we need to consider a quantity that remains relevant on 
both sides of the interface, and there are no chemical gauge 
issues concerning the momentum of the superfluid (essen- 
tially since it follows from the phase of the macroscopic 
quantum wavefunction). However, we still have to be care- 
ful. Basically, the presence of vortices and potential pinning 
complicates the picture. To make progress we will take the 
view that the irrotational condensate and the vortices can 
be considered separately. At the end of the day, the total 
momentum equation involves an average over these compo- 
nents. In effect, the interface must reflect this large scale 
average. Our approach to the problem represents this, yet 
it is admittedly rather naive. A number of issues need to 
be better understood, in particular concerning the way in 
which vortices extend from the fluid core to the elastic envi- 
ronment of the crust. Let us simply mention two problems: 
First of all, we know from low-temperature laboratory su- 
perfluids that vortices connect orthogonally to solid walls. 
It is natural to ask if the same is true for vortices that pen- 
etrate the neutron star crust-core interface. Secondly, the 
entrainment in the core leads to the vortices being magne- 
tised, due to the entrainment, but this effect relies on the 
protons being superconducting so is not active in the crust. 
The upshot is that a magnetised v ortex somehow connects 
to an unmagnetized one. How does this work? The answer 
may be linked to the transition from superconducting pro- 
tons to ones locked in nuclei. In this case, one would expect 
the presence of a current sheet. Presumably, this may also 
resolve any issues concerning the magnetic vortices, but the 
details have not yet been considered. These and other issues 
need to be resolved by future work. In the following we will 
adopt the pragmatic view that really difficult problems are 
perhaps best ignored. 

Anyway, considering first the irrotational part, we have 

dtPi + Vi (^Ax - ^v^ + vip]^ = (113) 

where x — [n, f] depending on whether wc are in the core or 
the crust. Working in the crust frame, and integrating over 
a small volume, this leads to the interface condition 

It should be noted that, in the particular choice of chemical 
gauge where we consider all the neutrons in the crust we 
have f — >^ n. In that case, the neutron chemical potential 
must be continuous, as expected. 

Moving on to the vortices, wo still do not have to worry 
about chemical gauge issues. Basically, one would expect 
each vortex to penetrate into the crust leading to the vortex 
density riv being continuous. Hence, we consider (13), where 
we recall that the vortex velocity depends on whether 
there is pinning or not. Working in the crust frame (as be- 
fore) and integrating, we arrive at the condition 

{N,W[)wr = iN,W^)w':^ (115) 

where we have used the fact that (in the present analysis) 
we do not consider pinning in the core. However, because 
of the chemical gauge invariance, we have Wf = Wi which 
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means that our final interface condition is 



(116) 



This result is extremely intuitive. If the crust vortices are 
free, moving with v|, then we must have v} = at the 
interface. Meanwhile, if the vortices are pinned, moving with 
vl, then the condition should be vl = v^. 

These results rojjrcsciit two limiting cases. In the general 
case, we would also need to keep track of the vortex velocity. 
This is, in principle, straightforward but we will not discuss 
the results here. It may be worth noting that such models 
may be considered from a three-fluid point-of-vicw, with the 
vortices forming a distinct "species". This is an interesting 
strategy that could prove advantageous in some situations. 



5.3 The perturbed problem 

Let us now move on to the conditions that need to be im- 
posed at the linear perturbation level. We can think of two 

possible strategics. Either wo take the view that the con- 
ditions derived in the previous section arc "exact" , which 
means that we can perturb them directly, or we start from 
the relevant perturbation equations and carry out the analy- 
sis for a small volume straddling the interface all over again. 
In principle, these two approaches should lead to the same 
answer. In practice, we find it useful to use a combination 
of them. 

We begin by considering particle conservation. In gen- 
eral, Lagrangian variation of the continuity equations leads 
to 

{dt + £„J (Axfix + nxVi^i) = (117) 
We require the trivial solution, i.e. take 



AxTix + n^Viii = 



(118) 



In the case of the protons, we can integrate this equation 
over the small cylinder across the interface. Taking the vol- 
ume small enough that variations in rip can be neglected 
(we are not allowing for density discontinuities), we then 
find that 



(119) 



As one would have expected, the normal component of the 
proton displacement should be continuous. This condition 
follows immediately if we want to avoid there being a void 
(or overlap) in the proton fiuid at the interface. 

To get the second condition , we perturb the equation 
for total baryon conservation. This leads to (in the crust) 

Acn = Ac(nf-hnc) = -ricV jm-V j {nt£,i^ (120) 

If we (again) assume that the densities are smooth, then we 
only need to consider 



Acn w -Vj (ucQ: -h nf^^^ 



(121) 



After integration across the interface, we find that we must 
have 

[nn (d - ^L) - nt (a - ) ] = ° (^^^^ 



In the case of the comprehensive gauge, when nt = rin, this 
reduces to 



AT, 







(123) 



Again, this condition seems natural, and we learn that the 
more complicated nature of (122) results from the fact that 
not all neutrons in the crust are free. We also see that we 
must keep careful track of the different number densities 
across the interface. 

Next we need the perturbed versions of (106) and (107). 
To derive these, we assume (as in the standard level set 
method) that the general conditions can meaningfully be 
extended away from the interface, and be perturbed in the 
usual way. There may be some technical issues associated 
with this approach, but we will not go into the details of 
this here. To derive the relevant conditions we start from 
(104) which leads to 



N. 



= Nj [p„(l - e^)w^^dtid - ^i) + SlA^p - ApT,^] . 

(124) 

It is easy to see that, for an axisymmetric background, the 
first term in each expression does not contribute to the nor- 
mal component. Contracting with TV' we are left with 

[AcP - N'NjA^T, ^1 = [ApP - N'NjA^T, ^1 . 

L J crust L J core 

(125) 

That is, we arrive at the expected traction condition. The 
horizontal result is (obviously) more complicated. A projec- 
tion orthogonal to leads to 

N, L(l - st)w^Me, - e) ~ (5l - iV'iVOAcT, ^1 

L J crust 



U(l - £.)wt^dt{a - ii) - {Si - N'Ni)A^Ti 

L J cor' 



(126) 

As in previous cases, the additional complications arise from 
the choice of chemical gauge. 

In the magnetic field case, we need (51) and (52). It is 
also useful to note that (110) leads to 

{NiAcB*} = . (127) 

Prom this we see that 

{N'AcBi) = {B'N^ [Vi^] + Vj^f) ) . (128) 
We then have 

NiAcp-NjAcTi' 

Ni (a„p+^A,B'']-^Nj (B'AcBi + SiAcB") . 



87r 



47r 



(129) 



The vertical component becomes 



{Acp + ^^cB') = — ( N,B' ) {N'AcBi) , (130) 



while the horizontal condition becomes 



_L 



B^ AcBi + BiA^B^ 







(131) 



At the perturbation level, we also need to consider the 
elastic problem. As before, we assume that the background 
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configuration is relsixed (and the core is fluid!), in which case 
we have 



and 



In this case, the vertical condition becomes 

(AcP+|/iV,e^H2(AiV»iV^ViCj) , 
while the horizontal one can be written 

±" (/iiv^'(v,C'+v,.c?)>=o. 



(132) 
(133) 
(134) 

(135) 



Combining the elastic and magnetic results to arrive at 
the conditions to impose in the general case is, of course, 
straight forward . 

Finally, we consider the conditions relating to the su- 
perfluid component. Perturbing the scalar condition (114) 
we see that we should require 



AcAf - w^AcEf + ( i - £f ] Acwic 



Ap/in - WnpApEn + ^ ~ 



(136) 



As far as the vorticity is concerned, it follows naturally from 
the discussion leading up to (116) that we should have 



_Lj (Acvi - Ap<) = 



(137) 



This condition is (obviously) satisfied if we have -LJ (^j — 
^n) = in the free vortex case and _L} {^i — ^i) = when 
perfect jjinning prevails. This ensures that there are no kinks 
in the vortices at the interface. 



6 CONCLUSIONS 

We have developed a Lagrangian perturbation framework 
relevant for the conditions that apply in a mature neu- 
tron star, accounting for the presence of superfluid compo- 
nents, the elastic crust and the magnetic field. The consid- 
ered physics impacts on a wide range of astrophysical phe- 
nomena, from pulsar glitches to magnctar seismology and 
various gravitational-wave emission mechanisms. Hence, the 
reported theoretical developments provide us with a solid 
foundation to consider exciting applications. 

There is also significant scope for future improvements 
of the theory. Most importantly, the effort should be ex- 
tended to relativistic gravity. This would open the door to 
truly quantitative considerations of realistic neutron star 
models, e.g. based on a modern supranuclear equation of 
state with composition and thermal gradients. Develop- 
ments in this direction are in progress. 

We also need to improve our understanding of the 
physics involved. Our discussion highlights the need to know 
a wider set of parameters, like the superfluid entrainment 
both in the star's core and in the elastic crust. The various 
interactions involving superfluid vortices, from pinning to 
mutual friction, also need to be understood. While we can 



continue to advance our understanding of the phenomenol- 
ogy of these very complex systems, we need to impose realis- 
tic constraints on our models. This requires, if not a precise 
knowledge of the involved parameters, some idea of what 
the permissible ranges may be. To achieve this goal we need 
a continued dialogue across different branches of physics, a 
challenging but ultimately rewarding exercise. 
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